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In this paper we show that a normal total number-of-particle fluctuation can be obtained con- 
sistently from the static thermodynamic relation and dynamic compressibility sum rule. In models 
using the broken U{1) gauge symmetry, in order to keep the consistency between statics and dynam- 
ics, it is important to identify the equilibrium state of the system with which the density response 
function is calculated, so that the condensate particle number A^'o, the number of thermal depletion 
particles A'', and the number of non-condensate particles Nnc can be unambiguously defined. We 
also show that the chemical potential determined from the Hugenholtz-Pines theorem should be 
consistent with that determined from the equilibrium equation of state. The A''*''"^ anomalous fluc- 
tuation of the number of non-condensate particles is an intrinsic feature of the broken U{1) gauge 
symmetry. However, this anomalous fluctuation does not imply the instability of the system. Using 
the random phase approximation, which preserves the (7(1) gauge symmetry, such an anomalous 
fluctuation of the number of non-condensate particles is completely absent 
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I. INTRODUCTION 

The number-of-particle fluctuation {SN'^) in an equi- 
librium system is a fundamental statistic problem since 
its scaling with the number of particles {SN'^) cx N"' re- 
lates to the stability of the system. The fluctuation is 
called normal if 7 = 1 and anomalous if 7 > 1. In 
the latter case, it implies the system is unstable, since 
the isothermal compressibility ^ 00 in the thermo- 
dynamic limit (see Eq. ^ below). For example, for a 
non-interacting uniform Bose system below the critical 
temperature, the fluctuation of the number of conden- 
sate particles {SNq) cx N"^, and that of the number of 
non-condensate particles (SN^^) cx N^^"^ in the grand 
canonical ensemble, while {SN§) — (SN^^) cx N^^^ in the 
canonical ensemble Q; all are anomalous since 7 > 1. 
However, for a trapped ideal Bose gas, the fluctuation 
of the number particles is normal, since the confinement 
effectively suppresses the thermal fluctuation (2j . 

Recently the number-of-particle fluctuation of inter- 
acting Bose-condensed systems has attracted much theo- 
retical attention, but whether or not it is anomalous still 
has not been resolved, since different methods predict 
different values of7|l,|l|l|10,llll|l3IIl|l3. Par- 
ticularly, even for the Bogoliubov approximation, both 
7 = 4/3 H H S S and 7 = 1 la 113 scaHng laws 
have been obtained. 

In order to see how these controversies arise, it is use- 
ful to examine the bases of these model calculations. 
Refs. ^ j^] use an energy functional of the total num- 
ber of particles N, the number of thermal excited parti- 
cles Nex, and the single-particle spectrum ej:. Using this 
energy functional, the fluctuations of the number of con- 
densate and non-condensate particles can be calculated 
using the partition function in either the grand canoni- 
cal ensemble, canonical ensemble, or microcanonical en- 
semble, and a 7 = 1 scaling law was obtained for both 
the condensate and non-condensate number-of-particles 



fluctuations. One important observation, which is essen- 
tial to obtain the 7 = 1 scaling law in this approach, 
is that phonon excitations have been excluded from the 
single-particle spectrum. The 7 = 1 scaling law for the 
condensate fluctuation is also obtained in Ref. 
which a single-condensate-mode Hamiltonian is used. In 
Refs. 0,0;IE0|) the spectrum obtained by the Bogoli- 
ubov approximation was used and a 7 = 4/3 scaling law 
was obtained for the non-condensate number-of-particle 
fluctuation. However using the compressibility sum rule 
(see Eq. a j = 1 scaling law was obtained in Refs. 
0, ^3 for the total number-of-particle fluctuation in the 
same Bogoliubov approximation. From this brief litera- 
ture survey, it is understood that the number-of-particle 
fluctuation in an interacting Bose-condensate system is 
highly model dependent, since the definition of the con- 
densate fraction, the number of non-condensate particles, 
and the energy spectrum are highly model de pen dent. 
But the contradictory results in Refs. 0, Q, H, Q| with 
that in Refs. within the same Bogoliubov approx- 

imation deserve further investigation. 

As can be shown, the total number-of-particle fiuctu- 
ations of any equilibrium system can be calculated from 
the static thermodynamic relation 



N 



N 



keT dN 



N 



= pksTKT, (1) 



TO 



where N is the particle number operator with expectation 
value N, ksT is the temperature, /i is chemical poten- 
tial, is the volume of the system, and p = N/fl is the 
number density. On the other hand, the total number- 
of-particle fluctuation can be also determined by the fol- 
lowing dynamic compressibility sum rule 
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where x„„(q';u;) is the density response function. The 
number-of-particle fluctuations calculated from these two 
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relations must be consistent in any approximation. How- 
ever, we have seen that the static result 7 = 4/3 obtained 
in Refs. 0, H, 1^ is not consistent with the dynamic 
result 7=1 obtained in Ref . ■ This leads to a con- 

tradictory conclusion about the stability of the system 
since, as argued by Yukalov ^.n anomalous fluctu- 
ation of the number of non-condensate particles would 
inevitably lead the system to be unstable while an inter- 
acting Bose-condensed system is stable. 

We have seen that the above inconsistency of the stat- 
ics with the dynamics in the Bogoliubov approxima- 
tion is related to the separation of the condensate and 
non-condensate components when the broken Bose C/(l) 
gauge symmetry is used. In this case, the Bose field op- 
erator ip is split as 



(3) 



where ^{r) is the Bogoliubov order parameter, and 
S'ip{f), usually called the non-condensate field operator, 
represents both the dynamic excitation and thermal de- 
pletion out of the condensate. This subtlety in S^l! shows 
that the condensate and non-condensate components are 
strongly correlated, and the condensate component can 
not be treated just as a static reservoir. Instead, the dy- 
namics aspect of ^' must be taken into account in calcu- 
lation of the number of condensate and non-condensate 
particles for the purpose of calculating the number-of- 
particle fluctuation from Eq. To resolve the incon- 
sistency of statics with dynamics and to treat 4" as a 
dynamic quantity, it is crucial to identify the equilibrium 
state with which the density response function Xnn is 
determined. Using this equilibrium state as a reference, 
the number of condensation particles A'o, the number 
of dynamically excited particles Nnc, and the number 
of thermally depleting particle N can be unambiguously 
defined. 

We will also show that the anomalous fluctuation of 
the number of non-condensate particles is an intrinsic fea- 
ture of the broken U{1) gauge symmetry. It is well known 
that a direct consequence of the above broken L''(l) gauge 
symmetry is that the poles of the single-particle Green 
function defined with Sip and the total density response 
function are identical. Therefore, the fluctuation of the 
number of non-condensate particles inevitably follows the 
iV*/^ anomalous law, since the momentum distribution 
of this non-condensate particle always has a. sin- 
gularity in the long- wavelength limit, regardless of at 
what level the interacting Bose Hamiltonian is truncated. 
However, we shall show that this anomalous fluctuation 
does not imply the instability of the interacting Bose sys- 
tem. This is because, as we shall show, the condensate 
and non-condensate components can not be treated as 
linearly independent terms so that an anomalous fluctu- 
ation of the number of non-condensate particles does not 
necessarily indicate that the system is unstable. 

One way to avoid such an anomalous fluctuation of the 
number of non-condensate particles is to work with an en- 
semble in which the gauge symmetry is not broken |l3j |. 



so that we can avoid the entangling of particle-conserving 
collective excitations and single-particle excited state in 
the pole of the non-condensate single-particle Green func- 
tion. Indeed, using the random phase approximation 
with inclusion of exchange (RPAE) developed by Min- 
guzzi and Tosi 14], which keeps the U{1) gauge symme- 
try, we are able to show that while the total number- 
of-particle fluctuation is normal and consistently deter- 
mined from statics and dynamics, the anomalous fluctu- 
ation of the number of non-condensate particles is com- 
pletely absent. 

This paper is organized as follows. In Sec^I briefly 
summarize the rules to build the non-condensate single- 
particle Green function and the density response func- 
tion with the broken C/(l) gauge symmetry. In Sec. IIIII 
we examine how the consistency between Eqs. ^ and 
(0) can be obtained in calculating the numbcr-of-particle 
fluctuation, and interpret the physical meaning of the 
anomalous fluctuation of the number of non-condensate 
particles in the Bogoliubov approximation. In Sec. IIVI 
we carry out a calculation in the random phase approxi- 
mation with inclusion exchange and in the dielectric for- 
malism to support our interpretations presented in Sec. 
IIIII The discussions and conclusion are presented in the 
last section. 



II. SINGLE-PARTICLE GREEN FUNCTION 
AND DENSITY RESPONSE FUNCTION WITH 
BROKEN t/(l) GAUGE SYMMETRY 

In this section, we briefly summarize the rules to con- 
struct the single-particle Green function and density re- 
sponse function with broken U(\) gauge symmetry. The 
details can be found in Refs. [Idi Ha. Il7|. 

We start from the Hamiltonian for a homogenous in- 
teracting Bose system in the second-quantized form 



^ k k k ' 9Q 
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where a contact two-body potential with strength g is 
fi — Ej: — fi is the non-interacting single- 
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particle energy with respect to the chemical potential, 
and at and are creation and annihilation operators 
for the interacting Bose particles. 

Now using Eq. (jS)) in its momentum space form, i.e., 
replacing aj and uq with \/Nq, where A'o is the number 
of particles condensed onto the ground state fc = 0, one 
obtains the approximated Hamiltonian 
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where p = Nq/VI is the condensate density, = at + a^, 
and 



E 



(6) 



is the density operator for the non-condensate particles. 
The total density operator is 



(7) 



Equation jSJl provides the starting point for many ap- 
proximations in which the broken U{\) gauge symmetry 
is used. 

The single-particle Green function matrix defined with 
(5V) is [Hill 



Gap{q 7^ 0,t) = -(T^a^a(T)at 
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where 
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Solving the Dyson equation involving a 2 x 2 matrix self- 
energy Ect/3, the sin g le- p article Green function Gap has 
the general form [ij, |lg 



Gap{k) 

where 
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{aiujn + £j:)Sa/3 + ctp^-a-pik) 



D{k) 



(10) 



{k)] [ 



IUJ„ 



+ S]+_(fc)E_+(fc). 



(11) 



Here notation fc = (fc; iwn) is used. Various truncations of 
the Hamiltonian (jSJ correspond to select certain types of 
self-energy diagrams in such a way that the Hugenholtz- 
Pines theorem 18] 
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is satisfied, so that the pole determined by 
gapless in the long- wavelength limit. 

The density response function defined as 

XnTi(g;T) = -(T^/5^(t)/5_^(0)) 

can be written as ^3 

Xnn{q) = ^ t^a{q)Gafi((i)^{i{(l) + 

Q/3 



(12) 
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where Aq, is the vertex function describing process of 
(de)excitations (in)out of the condensate and x™ is the 
regular part response function. It is not obvious in this 
approach to identify the equilibrium state with which 
the density response function %„„ is determined. We 
shall show in the next two sections that the identification 



of such an equilibrium state is important to unambigu- 
ously define the equilibrium condensate particle number 
A^o, thermal depletion particle number A^, and the cor- 
responding thermal depletion single-particle Green func- 
tion, which are used to build Aq,, Sq^, Ga/j, and x^, 
and to calculate the non-condensate particle number Nnc-, 
so that the consistency between statics and dynamics in 
calculating the number-of-particle fluctuation can be ob- 
tained. 



III. NUMBER-OF-PARTICLE FLUCTUATION 
IN BOGOLIUBOV APPROXIMATION 

W now reexamine the number-of-particle fluctuation 
problem in the Bogoliubov approximation at finite tem- 
perature. A comment is deserved: even though we work 
at finite temperature, there are no thermal depletion par- 
ticles. 

The vertex function is A^ = and the self-energies 
are 

S++(g;ia;„) = S — {q\iuJn) = 2gpo, (15) 
S+_(g; iujn) = S-+(<f; iuJn) = gpo, (16) 

and the chemical potential determined by Hugenholtz- 
Pines theorem is 

Ai = S]++(0;0)-E+_(0;0)=,gpo- (17) 

Substituting the above S^^, Aq, and /i into Eq. H10|l . one 
gets the corresponding single-particle Green functions for 
the non-condensate particles [1^, Il7| 
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and 
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Here A = — gpo has been introduced for future conve- 
nience, and A = for temperatures below T^. 

Substituting the vertex functions and the Green func- 
tion into Eq. I|14|) , one gets the density response function 
Xnuiq^^^n) for the interacting Bose gas 



XnniqJ^n) 



P0£q ( 1 



1 



(22) 



4 



We first calculate the number-of-particle fluctuation 
from dynamics. Substituting Eq. (|22|l into Eq. one 
gets 



{6N 
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where cb = ygPoT^ is the Bogoliubov phonon velocity. 
Therefore, we get a normal number-of-particle fluctua- 
tion from the dynamics. 

Now let's use Eq. to calculate the number-of- 
particle fluctuation. The total number of particles is 



given 

where Nnc is calculated as 
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Direct calculation shows that 
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The anomalous N'^^'^ behavior of this equation can be 



-I 



seen, since in the long- wavelength limit, ui 

and ~ i, so the integrand has a singularity. 

On the other hand, using the chemical potential /x 
gpo, one gets 



N dfi 
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The number-of-particle fluctuation is the sum of Eqs. 
(I26|l and H27() . which is clearly not consistent with Eq. 
This is the inconsistency of the result in Refs. 
S [13 with that in Refs. H El- We should re- 
mind ourselves that Eq. J^l is a thermodynamic relation 
for an equilibrium system which is described by a set of 
equations of state. Also, according to flnite-temperature 
linear-response theory, the density response function is 
calculated from an equilibrium state. Of course the equi- 
librium states that are used in Eqs. Q and ^ should be 
the same. So what is the equilibrium state for the Bogoli- 
ubov approximation at flnite temperature? To find the 
answer, we notice that Eq. H27() is identical to Eq. H23|) 
and we get it from relation /i — gpQ. Therefore, the equi- 
librium state in the Bogoliubov approximation is identi- 
fied to be a state that all particles occupy in the fc = 
level and its equation of state is given by p = gpo = gp. 
This identification is sound since the relation /i = gpo is 
exactly the time-independent Gross-Pitaevskii equation 
for a uniform system without thermal depletion parti- 
cles. Indeed, as proved by Leggett 0], the Bogoliubov 
Hamiltonian can be obtained from the Hamiltonian Q 



by keeping the terms that have non-zero expectation val- 
ues in a subclass of states built from Gross-Pitaevskii 
ground state. If we adopt this identification: 

(i) The number-of-particle fluctuation from statics is 
given by Eq. (|27|l . which is now completely consistent 
with Eq. (|23fl . Both are normal, and therefore the system 
is proved to be stable, as it should be. 

(ii) Nnc is the number of particles excited out of the 
condensate due to its oscillation, i.e. it is the depletion 
of the Gross-Pitaevskii equilibrium state A^o in Eq. 
(|^ should be replaced by Nq, 



Nl, = No- Nn 



(28) 



the number of particles remaining in the condensate after 
Nnc particles are dynamically excited out of the conden- 
sate. The single-particle Green function Gaf3 describes 
the dynamic process of the oscillation of the condensate. 
With the broken U(l) gauge symmetry, the oscillation of 
the condensate is interpreted as ejecting particles, which 
become the non-condensate particles. This interpreta- 
tion of the single-particle Green's function gives a reason- 
able explanation of the density response function given 
by Eq. H22|l . It also explains the physical meaning of the 
result given by Eq. H2t)|l . We argue that, the Nnc dynami- 
cally excited particles form a non-interacting system with 
chemical potential given by A = 0. Equation ^ then 
represents the number-of-particle fluctuation of this non- 
interacting system. In order to see this, we calculate the 
density response function of this non-interacting system 
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A similar result with Eq. 129() is obtained by Meier and Zwerger [3 by calculating the phase fluctuation of the order 
parameter ^(r) of Eq. It is important to point out the difference of the physical meanings between Eqs. (I22|l and 
(|29|l . According to Eq. (jSJ, the number-of-particle fluctuation of the non- interacting system is 

"£K + l)+'4«ij: (30) 
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where = {e^^k — 1) ^. Here lj^ = y^£j(e^ + 2,9^0) 
is the Bogoliubov mode. This is exactly the same as 
Eq. (7) in Ref. obtained by Gior gini et. al. for the 
fluctuation of non-condensate particles {6N^^)/N. We 
notice that the leading order terms of this result is iden- 
tical to those of Eq. H26|l in fc ~ region, where the 
anomalous behavior arises, since « nr. Therefore, 

for this non-interacting system, the number-of-particle 
fluctuations obtained from statics and dynamics are also 
consistent, even though they are anomalous. However, 
this anomalous fluctuation is not an implication of insta- 
bility of the interacting Bose gas, since we have proved 
from both statics and dynamics that the total number-of- 
particle fluctuation is normal. This can also be seen by 
substituting Eq. H24|l into Eq. but replacing iVo with 
Nq] the anomalous fluctuation due to Nnc is completely 
canceled out. This calculation clearly shows the impor- 
tance of the dynamic aspect of the condensate reservoir. 



(iii) There is a new consistency. The chemical potential 
as a functional of the total number of particles N and the 
equilibrium number of particles iVo in condensate can be 
determined both dynamically from the Hugenholtz-Pines 
theorem (I17|l and statically from the equilibrium equa- 
tion of state. These two must be consistent with each 
other. However, there is a deeper physical meaning of 
this consistency. The equilibrium state described by the 
equation of state has a definite number of particles (here 
is Nq). Therefore, this new consistency shows that the 
Hugenholtz-Pines theorem is to restore the conservation 
of the number of particles. Indeed, it is well known that 
the Hugenhotz-Pines theorem is required by the continu- 
ity equation 20, 21|. 



In the next section, we shall show that these inter- 
pretations about the number-of-particle fluctuation, N^, 
Nnc, and the single-particle Green function Ga/j, as well 
as the Hugenholtz-Pines theorem in the Bogoliubov ap- 
proximation remain valid at the level of approximation 
in which all the terms in Eq. |(SJ) are kept. As examples, 
we consider the random-phase approximation with in- 
clusion of exchange (RPAE) developed by Minguzzi and 
Tosi 01 and the dielectric formalism by Fliesser et. al. 
[2^ . We shall not give the detailed derivations, since 
they are available in the literature. The steps presented 
here highlight the physics at hand. 



IV. RANDOM-PHASE APPROXIMATION AND 
DIELECTRIC APPROACH WITH INCLUSION 
EXCHANGE 

In the RPAE, the equilibrium equations of state of the 
Bose-condensed system are the time-independent flnite- 
temperature Gross-Pitaevskii equation for the conden- 
sate and static Hartree-Fock equation for the thermal 
depletion particles. For a homogenous system, they are 

= gpo + "^gp, (31) 

hHF^j:ir)^sS^ibf:{r), ^ ^ 0, (32) 

where 

hHF = -T, \-2gp~fi, (33) 

2 m 

=ej^ + 2gp-^i, (34) 

are the static Hartree-Fock Hamiltonian and single- 
particle energy with respect to the chemical potential 
for the thermal depletion particles. Here po = Nq/H, 
p ~ N/n, and p = N/D, are the equilibrium condensate, 
thermal depletion, and total densities with A'o, N and N 
the corresponding equilibrium condensate, thermal de- 
pletion, and total number of particles. We emphasize 
again that the number of particles in this equilibrium 
system is conserved. By neglecting the thermal deple- 
tion N, we arrive at the equilibrium equation of state for 
the Bogoliubov approximation. We notice that single- 
particle orbits for the condensate and thermal depletion 
particles are governed by two different Hamiltonians, and 
therefore, they are not generally orthogonal. However, 
for a uniform system, these single-particle orbits are sim- 
ple orthogonal plane waves. We also notice that there is 
a gap in the single-particle spectrum 

limef^=gpo. (35) 

This gap has important effects on many properties of a 
Bose-condensed system . We will come back this issue 
in the last section. 

In this paper, we concern its effect on the stability of 
the system, as we shall show in the following. 

We deflne a thermal depletion single-particle Green 
function for the static Hhf 

GHF{k;iuj,i) = (36) 

tuJn - 
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The number of thermal depletion particles is found to be 



N : 



1 .r-^ ~ ^ ft 
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where At = yj2^Flmk^ and z = e'^^^^^f'') = e-l^^po 
and g'y{z) is the Bose function. Equation H37I) is the 
equation of state equivalent to Eq. H32I) for the thermal 
depletion particles. The self-consistent relations among 
N , Nq, N and fi are given by Eq. H31|l and 

7V(/i) = Noifi) + Nifi) = 7Vo(m) + ^53/2[^(m)]- (38) 

At" 



itself is identically zero. Similar calculations for the RPA 
without the exchange show that (SN^) follows the A^^/"^ 
anomalous scaling law, but the total number-of-particle 
fluctuation is 



{SN 



(45) 



which is normal. 

We now calculate the density response function around 
the above equilibrium state and calculate the number-of- 
particle fluctuation from dynamics. 

The linearized equations for the density fluctuation 
have the matrix form [iJl 



The number-of-particle fluctuation can be calculated 
by substituting Eq. into Eq. Using the equations 
of states (|ST|) and l|?7|) . we find 



dNo n dN 
-1^^ + 2—, 

dN (3 , , fr. 



and as a consequence 
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where Pn is defined as 
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Equation (|40|l reduces to Eq. (|27() , obtained in Bogolibov 
approximation when one sets Pq = and po — p. 

We must also point out that ksT^ ^ {SN§), 

ksT^ ^ {6N^). Because of the ensemble used in the 



RPAE |2g, (SN^) and (SN^) are actually 
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These two results are clearly not the same those given by 
Eq. Therefore, 



{5N^) ^ (SN^) + (SN^). 



(44) 



This shows that, even in the mean-field level approxi- 
mation, the condensate and the thermal depletion com- 
ponents are strongly correlated. This is not surpris- 
ing, since below the critical temperature, the presence 
of the condensate pins down the chemical potential to 
be fi — 2gp + gpo and Eq. is a self-consistent rela- 
tion between N and p. Even in the Bogoliubov approx- 
imation, it is this self-consistent relation that predicts a 
number-of-particle fluctuation given by Eq. (|27|l consis- 
tent with Eq. H23|) while the fluctuation of the condensate 



6p 



Xcc "Kch 
Xhc X.hh 



SU" 



(46) 



where SW^ and SU" are the spatially and time-varying 
external potentials for the condensate and thermal de- 
pletion components. We emphasize here that the matrix 
form H46() of the density response function is simply a 
result by splitting the system into a condensate and a 
thermal depletion component, in which the number of 
particles is conserved, instead of a result of using the 
broken 11(1) gauge symmetry Q as claimed by Minguzzi 
and Tosi . The total density response function is then 
given by 



Xn 



Xcc + Xcfi + Xnc + Xf. 



(47) 



On the other hand, according to the linear response 
theory, 

5p^ = xl5Ujjp - xl{5U^ + gSpo + 2g5p), 
SP = xlUfjp = x^('5t/" + 2g6po + 2gSp), (48) 

where Xc ^^'^ Xfi the density response functions of the 
condensate and the thermal depletion around the equilib- 
rium state, respectively. From the above two equations, 
the four components are found 

Xcc - (1 - 2gxl)D-\l, xnn = (1 - gxl)D-\l, 
Xcn = 2gxlD^^x%, Xnc 
where 

D = {\- axl)(l - 2.gx°) - ^g\lxl (50) 

For homogenous systems, the density response func- 
tions of the condensate and of the thermal depletion com- 
ponent can be obtained by linearizing the time-dependent 
Gross-Pitaevskii equation for the condensate around Eq. 
(|31|l . and the time-dependent Hartree-Fock equation for 
the non-condensate around Eq. (|32|l . respectively. They 
are given by 



2gxlD-\l, (49) 



X°c{?\^^-n 
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where = (z~^e'^^£ — 1)"-^ is the occupation number spouse function for the homogenous Bose-condensed sys- 
of the static Hartree-Fock single-particle level of the ther- tem is 
mal depletion particles. Therefore, the total density re- 
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xliq; iujn) + 2pe,- [i + gxliq; i^n)] 


iiujn)^ - el 


[1- 





(53) 



Now substituting Eq. 153|) into Eq. |2Jl, one gets the 
total number-of-particle fluctuation from dynamics 



N p mc% 1 + 2gPo ' 



(54) 



of the fact that 



Here we have made use 

lim|g-|^oXs('7'.0) = -^9i/2{e-^^^^) = n as given 
by Eq. (|^ . One can see that Eq. is exactly the 
same as Eq. (gU]). 

We have shown the total number-of-particle fluctua- 
tion is normal and consistent between statics and dy- 
namics in the RPAE, and therefore, the interacting Bose 
system is proved to be stable. In the above derivations, 
the number of particles is conserved and there is not any 
anomalous number-of-particle fluctuation. This is be- 
cause in this RPAE, the numbers of particles in the con- 
densate and thermal depletion component are not time- 



dependent jl9i | and do not change with the external po- 
tential because of entropy conservation Therefore, 
they always take the equilibrium values. The above re- 
sults can be derived in a more general time-dependent 
Hartree-Fock scheme which preserves the U{1) gauge 

symmetry r23:i. 

Now in order to see how the anomalous fluctuation 
arises when the broken £7(1) gauge symmetry is used, we 
can follow the steps in Ref. [221 to build the vertex func- 
tion Aq,, self-energy Eq,^, and the single-particle Green 
function Gap by using using Eqs. and as the 
reference. This means that the equilibrium condensate 
A'o, thermal depletion TV, and Eq. (|36|l should be used to 
build up Aq, Sq^, and regular x^, not the as yet to be 
determined iVg, Gaf3, and Nnc- Here we skip those steps 
and just cite the final results for the single-particle Green 
function matrix below 



J 



G++{q;iuJn) = G — {q;-iujn) = 



G+-iq;iuJn) = G-+{q;iuJn) 



[iun - eg] [l - 2gxl{q; iuJn)] + gpQ [l + gxUq; "i^n)] 



el 
g 



[1 - 2gxfM *^«)] - 2poeg [1 + 2gx°M «^")] ' 
gpo [l -hgx°(9;iwn)] 



(i^„)2 - el [1 - 2gxl{q;iLo^)] - 2poeg [1 + 2gxl{q;i^n)] 



(55) 



and the density response function x„„((f; iwn), which is 
same as Eq. (|^ . 

Since both the equilibrium state and the density re- 
sponse function are the same as in RPAE, therefore, one 
gets the same consistent normal number-of-particle fluc- 
tuations from statics and dynamics in this dielectric for- 
malism as those in the RPAE j23| . 

The chemical potential from the Hugenholtz-Pines the- 
orem in this approximation turns out to be 

PHP = S++ (0; 0) - !]+_ (0; 0) = gpo + 2gp, (56) 

which is the exactly the same as Eq. (|31|l . 

Using G-i--i-{q;iu!n), the number non-condensate par- 
ticles Nnc and its fluctuation {6N^^) are found to be 



the same as Eqs. (1^ and l|Sn|) . For a dilute gas 

pa^ <C 1, where a — P is usually very small because 
of the single-particle gap [2^. It is thus straightforward 
to show that the single-particle Green function Ga/s given 
by Eq. H55|) has the similar form as that in Bogoliubov 
approximation for small k. For example, the pole is given 



by ujr fa csk 



1 



Therefore, foUow- 



2gP{k,uj^ csk) 

ing the steps as in the Bogoliubov approximation, one 
can show that {6 Nnc) follows the 7 = 4/3 scaling law. 

Since N from Eq. H37|) is the number of thermal deple- 
tion particles, the difference 



5No = Nn 



N 



(57) 
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can be interpreted as the number of particles excited out 
of the condensate by the oscillation of the whole system 
induced by the external time-dependent potential. In- 
deed, in the Bogoliubov approximation SNq — Nnc since 
the depletion of the condensate is completely caused by 
the dynamic collective excitation. Therefore, the single- 
particle Green functions H55|) can be interpreted as dy- 
namic ones comparing to the thermal depletion Gh f de- 
fined by Eq. H36|l . As in Bogoliubov approximation, this 
interpretation is allowed only because of the broken C/(l) 
gauge symmetry. The total number of particles is ex- 
pressed as 

N = N'^ + 7V„„ (58) 

where Nq = Nq ~ SNq. When substituting Eq. (|5H|l into 
Eq. to calculate the number-of-particle fluctuation, 
the anomalous fluctuation due to Nnc is exactly canceled 
out by the one from TVq, so that the total number-of- 
particle fluctuation is normal, which is given by Eq. (|40|l . 

We can see that the anomalous fluctuation of the num- 
ber of non-condensate particles {SN^^) is solely due to the 
single-particle Green functions defined by Eq. I|10() whose 
poles entangle the single-particle and particle-conserving 
collective excitations, a directly consequence of the U{1) 
symmetry breaking rather than an implication of the in- 
stability of the Bose system since the total number-of- 
particle fluctuation is normal and consistent from statics 
and dynamics. More than thirty years ago, Straley ad- 
vised caution psj in using such a single-particle Green 
function to describe the zero-sound characteristic spec- 
trum of the superfluid ^He. Leggett also argued that 
there are no circumstances in which Eq. ||2Jl is physically 
correct. 



V. DISCUSSIONS AND CONCLUSION 

We have shown that the anomalous fluctuation of the 
number of non-condensate particles is an intrinsic feature 
of the broken U{1) gauge symmetry and is completely 
absent in the RPAE in which the U{1) gauge symmetry 
is preserved. This may be just related to the different 
interpretations of the dynamic process of the condensate 
of these models, as we point out in previous sections. But 
since this anomalous fluctuation of the number of non- 
condensate particles has not any physical significance, 
we can safely say that it is just a by-product of using the 
broken U{1) gauge symmetry. 

An advantage to keep the C/(l) gauge symmetry is that 
the single-particle spectrum and the collective excitation 
spectrum are distinguished. As shown in the RPAE, they 
are the poles of thermal depletion single-particle Green 
function and the density response function, respectively. 
But single-particle spectrum has important effects on the 



collective mode. For example, because of the gap, Xh is 
quite small for a dilute gas so that the thermal depletion 
component does not sustain a zero sound mode. The 
whole effect of the thermal depletion particles is to shift 
the Bogoliubov mode. Also the gap causes the damping 
of the thermal depletion particles to the collective mode 
exponentially decreases when the temperature drops. As 
we recently discussed 631, this single-particle gap is also 
responsible for the completely screening of the external 
potential by the condensate. 

By using the broken U{1) gauge symmetry, the bound- 
ary of the single-particle spectrum and the collective 
mode are entangled together if the poles of Ga/s are in- 
terpreted as the single-particle excitations. But as far 
as the dielectric formalism in Ref. concerned, there 
is a gapped single-particle spectrum as same as that in 
the RPAE for the equilibrium reference. In this sense, a 
gapped single-particle spectrum and a gapless collective 
mode do coexist even in the dielectric formalism. This 
is another point to identify the equilibrium state with 
which the density response function is calculated. 

As pointed out by Meier and Zwerger Jj, the anoma- 
lous fluctuation of the non-condensate particles is related 
to the gapless mode in the superfluid Bose system. Our 
analysis shows this is the case if one uses the broken Bose 
U{1) gauge symmetry. In the RPAE, which preserves this 
gauge symmetry, there is no such an anomalous number- 
of-particle fluctuation related to the gapless superfluid 
mode. 

In conclusion, we have shown that in models using 
the broken C/(l) gauge symmetry, the number-of-particle 
fluctuation is normal and can be calculated consistently 
from the static thermodynamic relation and dynamic 
compressibility sum rule if the equilibrium states are 
identified. We also show that the chemical potential de- 
termined from the Hugenholtz-Pines theorem should also 
be consistent with that determined from the equilibrium 
equation of state. The TV''/^ anomalous fluctuation of 
the number of non-condensate particles is an intrinsic 
feature of the broken U{1) gauge symmetry. However, 
this anomalous fluctuation does not imply the instabil- 
ity of the system. Using the RPAE, which preserves the 
U{1) gauge symmetry, such an anomalous fluctuation of 
the number of non-condensate particles is completely ab- 
sent. 
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